Next-to-leading order QCD corrections to the top quark decay 
via model-independent FCNC couplings 
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DO and CDF collaborations at the Fermilab Tevatron have searched for non-standard-model single 
top-quark production and set upper limits on the anomalous top quark flavor-changing neutral 
current (FCNC) couplings fi^/A and k?„/A using the measurement of total cross section calculated 
at the next-to-leading order (NLO) in QCD. In this Letter, we report on the effect of anomalous 
FCNC couplings to various decay branching ratios of the top quark, calculated at the NLO. This 
result is not only mandatory for a consistent treatment of both the top quark production and 
decay via FCNC couplings by DO and CDF at the Tevatron but is also important for the study of 
ATLAS and CMS sensitivity to these anomalous couplings at the CERN LHC. We find that the 
NLO corrections to the partial decay widths of the three decay channels t — > q + g, t ^ q + 7 and 
t — » q + Z are at the order of 10% in magnitude and modify their branching ratios by about 20%, 
0.4% and 2%, respectively, as compared to their leading order predictions. 
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PACS numbers: 12.15.Mm,12.38.Bx,14.65.Ha 

The top quark FCNC processes t -> q + V (V = 
g, 7, Z) have tiny branching ratios in the standard model 
(SM) and are probably unmeasurable at the CERN Large 
Hadron Collider (LHC) and future colliders. On the 
other hand, any positive signal of these rare decay events 
would definitely imply some new physics beyond the SM. 
As the LHC will produce abundant top quark events 
(about 10 8 per year), even in the initial low luminosity 
run (~ 10 fb~ /year) 8 x 10 6 top quark pairs and 3 x 10 6 
single top quark will be produced yearly, one may antic- 
ipate to discover first hint of new physics by observing 
anomalous couplings in the top quark sector. 

Recently, from their measurements of the total cross 
sections both DO [l[ and CDF collaborations at the 
Fermilab Tevatron have searched for non-standard-model 
single top-quark production and set upper limits on the 
anomalous FCNC couplings k? c /A and K,f u /A, where 
the leading order (LO) cross sections have been scaled 
to NLO 0] (or resummation 0]) predictions. At the 
LHC, ATLAS collaboration has presented its sensitiv- 
ity to studying FCNC top decays [5| . These results show 
that studying top quark FCNC couplings will provide 
a good probe to new physics beyond the SM. Although 
there are many discussions in the literature on rare decay 
processes involving model-independent top quark FCNC 
couplings d, 0, @, d, E3, U > 13 , most of them were based 
on LO calculations. However, especially for t — ► q + g, 
due to the large uncertainties from the renormalization 
scale dependence in its LO prediction through the strong 
coupling constant a s , it is necessary to improve the theo- 
retical prediction to NLO in order to match the expected 
experimental accuracy at the LHC. 

In this Letter, we present the analytic results of the 



NLO QCD corrections to the partial decay widths and 
decay branching ratios of top quark via anomalous FCNC 
couplings for the three processes t — * q + g, t — > g + 7 
and t — > q + Z, which are the most commonly studied 
decay channels by the experimentalists at the Tevatron 
and the LHC. 

New physics effects involved in top quark FCNC 
processes can be incorporated in a model-independent 
way into an effective Lagrangian which includes the 
dimension-5 operators as listed below [l(| 
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where A is the new physics scale, 9w is the Weinberg 



angle, and T a are the Gell-Mann matrices. 



^9 



with 



V = p,7 and Z are normalized to be real and positive, 
while fY q and \i( q are complex numbers in general. Since 
the partial widths in our calculation depend only on the 
product of k and |/| 2 + \h\ 2 , we could set |/| 2 + \h\' 2 = 1 
by redefining the parameter k. 

From the effective Lagrangian given by Eq. (Q]), we 
obtain the following LO partial decay width of the FCNC 
top decays in D = 4 — 2e dimension, 
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where the masses of light quarks q (q = u, c) have 
been neglected, flz = (1 — A/f/m 2 ) 1 / 2 and C e = 

r(2~2e) {^if) e - Here ' we also define r °(* Q + V) = 

r"o (i - * g + V")lf->n; which are consistent with the expres- 
sions in Ref. flOj . 

Below, we present in details our calculation for the 
inclusive decay width of the top quark, decaying into 
hadrons, up to the NLO with the LO partonic process 
denoted as t — ► q + g. The final results of t — > q + 7 and 
t — ► q + Z are also given for completeness. 

At the NLO, we need to include both one-loop virtual 
gluon corrections and real gluon emission contribution. 
We use D = 4 — 2e dimensional regularization to reg- 
ulate both ultraviolet (UV) and infrared (IR) (soft and 
collincar) divergences so that all divergences appear as 
l/e a with a = 1 and 2. The UV singularities cancel after 
summing up the contributions from the one-loop virtual 
diagrams and counterterms according to the same con- 
vention used in Ref. The soft-singularities cancel af- 
ter adding up the virtual and real radiative corrections. 
To cancel collinear singularities, we need to also include 
contribution induced from gluon splitting to light quark 
pairs at the same order in QCD coupling. 

The renormalized virtual corrections to the partial de- 
cay width of t — > q + g is 
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Here fi is a renormalization scale and Nf denotes the 
number of active flavors at the scale \i. All the ultraviolet 
divergences have canceled in r^ irt , as they must, but the 
infrared divergent piece is still present. 

The contribution from real gluon emission (t — > q + g + 
g) is denoted as T rca \(t — > q + g + g). In order to cancel all 
the collinear singularities in the sum of virtual and real 
radiative corrections, we also need to include the con- 
tributions from gluon splitting into a pair of quark and 
anti-quark in the collinear region, which is denoted as 
r re ai(i —> q + q' + q')- Note that there are two configu- 
rations of final state when the flavor of the light quark 
coming from gluon splitting is the same as the light quark 



directly from the FCNC coupling, and only one configu- 
ration when they are different. We find that at the NLO 
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r rC al(i -» q + q' + q') = -» q + g) 

07T 



—N f 



eiR 



(5) 



After adding them together, the total real contributions 
can be written as 
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Combining the real and virtual contributions, we obtain 
the full NLO corrections for t — > q + g as 



r N Lo(i-9 + .9) = H irt + r? eal 



72tt 



174 In 



"'7 



-127V/ In ( -^2 ) - 36/V/ - 38tt 2 + 749 



(7) 



Note that all the infrared divergences have been canceled 
in Eq. ([7]) as they should. 

As for t -»■ q + 7 and t -> q + Z, the NLO QCD 
corrections to the partial decay widths are given by 
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TABLE I: Branching ratios as functions of Kt q /A. Here \x 
m t . 



BR [in unit of (^TeV) 2 ] 


LO 


NLO 


NLO/LO 


t -> q + g 


1.0010 


1.1964 


1.195 




0.0544 


0.0542 


0.996 




0.0448 


0.0458 


1.022 
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respectively. Hence, up to the NLO, the partial decay 
widths of the three FCNC decays can be obtained by 

r(t -» 9 + v) = r (< -» q + v) + V N lo(* -> « + v)- 

In order to study the effect of NLO corrections to decay 
branching ratios, we define the following branching ratio 
for latter numerical analysis: 
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The partial decay width for the dominant top quark de- 
cay mode of t — > VF + at the NLO can be found in 
Ref. [l3j], namely, 
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where /3w = (1 - m^/m 2 ) 1 / 2 . 

For the numerical calculation of the branching ratios, 
we take the top quark width given in Ref. [l3| with the 



TABLE II: Decay widths as functions of kJ^/A. Here n — nit- 



Width [in unit of (^TeV) 2 GeV] 


LO 


NLO 


NLO/LO 


t -> q + g 


1.443 


1.577 


1.09 


t -> 5 + 7 


0.078 


0.071 


0.91 




0.065 


0.060 


0.93 
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FIG. 1: Branching ratio of t — > q + g as function of K B tq j 'A. 
The DO upper limits from Ref. [l[ is shown. 

following parameters 

(m w , m z , m t ) = (80.398, 91.1876, 171.2) GcV, 
Nf = 5, a = 1/128, sin 2 W = 0.231, 

V tb = l, Gf = 1.16637 x 10~ 5 GcV -2 . 

We analyze our results as functions of FCNC couplings 
kJ^/A and rcnormalization scale fi. For /i = m t , our 
numerical results are listed in Table U and Table |TT] show- 
ing the NLO effect to various decay branching ratios and 
partial decay widths, respectively. From Table QJ we see 
that the NLO correction increases the LO branching ratio 
by about 20% for t — ► q + g, while the NLO corrections 
are much smaller for the other two decay modes. But 
we note that the NLO results modify the LO widths by 
about 10% in magnitude for all the three modes. 

For convenience, we show the branching ratio of t — > 
q+g as function of n 9 tq / A in Fig.[TJ where we have set (i = 
m t as before. Using the upper limits measured by the DO 
Collaboration at the Tevatron [![, we get the following 
results: 



A 



< 0.15 TeV" 



< 0.037 TcV" 

A 



BR(t -^c + g) < 2.69 x 10~ 2 , 
BR(t -^u + g) < 1.64 x 10~ 3 . 



Following the analysis in Ref. [f| , we plot the coupling 
/A as a function of the branching ratio in Fig. [5] where 



1 tq 



the ATLAS sensitivities for the two different expected in- 
tegrated luminosities are also exhibited. From Fig. [5J it 
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FIG. 2: Coupling k^/A as function of branching ratio. The 
ATLAS sensitivities from Ref. [f| is shown. 
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the decay branching ratios and partial decay widths of 
the top quark. We define Rlo(^) = r o(/i) /I o(/x = rn t ) 
and i?NLo(/ i ) = r(/i)/r(/x = m t ), and show the value of 
R(/J.) as functions of /x for i — ► q + p in Fig. [3] It shows 
that the theoretical uncertainty from the renormalization 
scale dependence can be largely reduced to a couple of 
percent once the NLO calculation is taken into account. 

In conclusion, in order to achieve consistent studies for 
both the top quark production and decay via FCNC cou- 
plings, we have calculated the NLO QCD corrections to 
the three decay modes of the top quark induced by model- 
independent FCNC couplings of dimension 5 operators. 
For t — ► q + g, the NLO results increase the experimen- 
tal sensitivity to the anomalous couplings. Our results 
show that the NLO QCD corrections enhance the LO 
branching ratio about 20%. Moreover, the NLO QCD 
corrections vastly reduce the dependence on the renor- 
malization scale, which leads to increased confidence in 
our theoretical predictions based on these results. For 
t — > q + 7 and t — ► q + Z, the NLO corrections are mi- 
nuscule in branching ratios, albeit they can decrease the 
LO widths by about 9% and 7%, respectively 
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FIG. 3: Ratio 7? as function of the renormalization scale fi. 



is palpable that the NLO prediction improves the sen- 
sitivities of the LHC experiments to measuring the top 
quark FCNC couplings. With an integrated luminosity 
of 10fb _1 , the ATLAS experiment sensitivities can be 
translated to the following relations on FCNC couplings: 

K 9 

BR(t -> q + g) > 1.3 x 10" 3 => > 0.033 TeV -1 , 

7 

BR(t -> q + j) > 4.1 x 10~ 5 > 0.028 TcV~\ 

K Z 

BR(t -> q + Z) > 3.1 x 10~ 4 => > 0.082 TeV" 1 , 

and with an integrated luminosity of lOOfb -1 , they can 
be translated to the more stringent relations: 

K 9 

BR(t -► q + g) > 4.2 x 10~ 4 => > 0.019 TeV" 1 , 

7 

BR(t -> q + j) > 1.2 x 10" 5 -jt > 0.015 TeV" 1 , 

K z 

BR(t^q + Z) >6.1 x 10~ 5 => > 0.036 TeV" 1 . 

Lastly, we illustrate the fact that the NLO prediction 
reduces the theoretical uncertainty in its prediction on 
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